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New theory is presented for evolving systems, which are autonomously controlled 
subsystems that self-assemble into a new evolved system with a higher purpose. Evolving 
systems of aerospace structures often require additional control when assembling to 
maintain stability during the entire evolution process. This is the concept of Adaptive Key 
Component Control that operates through one specific component to maintain stability 
during the evolution. In addition, this control must often overcome persistent disturbances 
that occur while the evolution is in progress. Theoretical results will be presented for 
Adaptive Key Component control for persistent disturbance rejection. An illustrative 
example will demonstrate the Adaptive Key Component controller on a system composed of 
rigid body and flexible body modes. 


I. Introduction 

E VOLVING systems are autonomously controlled subsystems which self-assemble into a new evolved system 
with a higher purpose 1-2 . Evolving systems of aerospace structures often require additional control when 
assembling to maintain stability during the entire evolution process 3 * " 5 . An adaptive key component controller 
has been shown to restore stability in evolving systems that would otherwise lose stability during evolution 6 " 8 . The 
adaptive key component controller uses a direct adaptation control law to restore stability to the evolving system 
through a subset of the input and output ports on one key component of the evolving system. The reader is directed 
to refs. 7-8 for a more detailed description of evolving systems. In this paper, we address the situation where 
persistent disturbances can appear in some components and must be mitigated by the adaptive key component 
controller. Such disturbances will often be attendant in actively controlled rendezvous and docking or component 
capture applications. 

The control laws used by the adaptive key component controller to restore stability in an evolving system are 
guaranteed to have bounded gains and asymptotic tracking if the evolved system is almost strictly dissipative. 
Hence, it is desirable to know when the dissipativity traits of the subsystem components, including the key 
component, are inherited in an evolving system. We present results describing when an evolving system will inherit 
the almost strict passivity traits of its subsystem components. Then we will present an adaptive key component 
controller that restores asymptotic stability with bounded adaptive gains and mitigates the effect of persistent 
disturbances during evolution. The paper is organized as follows: the next section describes the framework for 
evolving systems, then the conditions under which dissipativity traits are inherited in evolving systems will be 
presented. Next, the design of adaptive key component disturbance accommodating control will be described. 
Finally, an illustrative example is presented to demonstrate the concepts presented earlier. 
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II. Evolving Systems Framework 


A. Mathematical Formulation of Evolving Systems 

A mathematical formulation of a nonlinear time -invariant evolving system is given here. Consider a system of L 
components of individually, actively controlled subsystems which can be described by the following equations for 
the ith component: 


\ X i ~ f i ( X i ’ U i ) (1) 

1 y i = g,(x i ,u i ) 

where / = 1. 2 L . The /' h component has a Performance Cost Function ./, and a Lyapunov Function Vj. These are 

the building blocks of the evolving System. When these individual components are joined to form an evolved 
system, the new entity becomes: 


fx = f(x,u ) 
\y = g(x,u) 


with .V = | .V| ] 7 . y = [_y r .._y L ] r , Performance Cost Function J, and Lyapunov Function V. The ; lh component 
in the above evolved system is given by: 


L 

*i = ft ( x i . ) + Z (x i ’ x j ’ u j ) (3) 

7=1 

where 0 < 8- < 1 and f ij (x i ,X j ,U j ) represents the interconnections between the i th and / h components. Note 
that when £- — 0 , the system is in component form and when £- — 1 , the system is fully evolved. As the system 
evolves, or joins together, the 8- ’s evolve from 0 to 1. 

The components of the evolving system are actively controlled by means of local control. Local control means 
dependence only on local state or local output information, i.e., U i = /z. (x. ) or U i = h ( (y.) . In general, the local 
controller on the / th component would have the form: 

U- =h(yi’Zi) 

where Z- is the dynamical part of the control law. Local control will be used to keep the components stable and to 

meet the individual component performance requirements, J ] . 

Once the system is fully evolved, the / th component in the fully evolved system becomes: 

L 

x i = fi ( x i V M, ) + Z fij {x ‘ ’ x j ’ U J ) (5) 

7=1 

A state space version of the / th individual component of an evolving system where the components are connected 
through the states can be represented as: 
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( 6 ) 


L 

A = A ( x , ) + B , ( x , X + Z % \ (x , . x j ); *,■ (0) = ^o, 

* ;=l 

ji=CA x i) 

where / = 1, 2, . . . L , X. = [x[ . . . X^ ] T is the component state vector, = [wj . . . ] T is the control input vector, 
v, =[v;-- •y p ] T is the sensor output vector, (A.(x.), 5. (x.), C-(x.)) are vector fields of dimension 
n { x n { , n. X m ? . , and p t X ft. , respectively, and the connection forces between components are represented in the 

ft. X ft ■ connection matrix, A..(x.,X.) with S ■■ — £■■ . The state space representation of the evolved system then 

1 j y * j x ji y 

becomes: 


fx = A(x) + B(x)u 

{ (J) 

b = c(x) 

which will also be written as (A(x), 5(x), C(x)) . 

B. Characteristics and control of evolving systems 

We say a subsystem trait, such as stability, is inherited when the evolved system retains the characteristic of the 
trait from the subsystem. Previous papers have examined the inheritance of stability and shown that stability is not a 
generally inherited trait in evolving systems 3-5 . Control can be used in evolving systems to maintain stability in a 
system that would fail to inherit stability during evolution. 

Key component control has been proposed to restore stability in evolving systems that would otherwise lose 
stability during evolution. The approach used for key component control design is for the control and sensing of the 
components to remain local and unaltered except in the case of one key component which has additional local 

control added to stabilize the system 
during evolution. The key component 
controller operates solely through a 
single set of input-output ports on the 
key component, see Figure 1 . 

Only the key component of the 
evolving system needs modification to 
restore the inheritance of stability. A 
clear advantage of the key component 
design is that components can be reused 
in many different configurations of 
evolving systems without the need for 
component redesign. The reuse of 
components which are space-qualified, 
or at least previously designed and unit 
tested, could reduce the overall system 
development and testing time and should result in a higher quality system with potentially significant cost savings 
and risk mitigation. 

In many aerospace environments and applications, the parameters of a system are poorly known and difficult to 
obtain. Adaptive key component controllers, which make use of a direct adaptation control law, are a good design 
choice for restoring stability in evolving systems where access to precisely known parametric values is limited. The 
sufficient condition for a nonlinear evolving system with an adaptive key component controller to be guaranteed to 
have bounded gains and to have asymptotic output tracking is that the system be almost strictly dissipative. So, we 
are interested in the conditions under which the inheritance of almost strict dissipativity can be guaranteed in 
evolving systems. The next section gives definitions and theorems related to dissipativity. 



Fig. 1. Evolving system with key component controller. 
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III. Inheritance of Almost Strict Dissipativity in Evolving Systems 

Inheritance of almost strict dissipativity of subsystems is desirable in evolving systems that use an adaptive key 
component controller to restore stability. Consider a nonlinear system given by Eq. 7. We say this system is Strictly 
Dissipative when there exists V (x) > 0 \/x ^ 0 such that the Lie derivatives satisfy: 

f L a V = VVA(x) < -S(x) \/x 
\l b V = VVB(x) = C t (x) 


where W = gradientV . 

The function V (x(t)) is called the Storage Function for Eq. 7, and the above says that the storage rate is always 
less than the external power. This can be seen by taking the time derivative of V (x(t )) to obtain: 

V = VV[A(x) + B(x)u] 

< -S(x) + C T (x)u (9) 

= -S(x) + (y,u) 


Taking u = 0 , it is easy to see that Eq. 9 implies the first equation in Eq. 8 but not necessarily the second equation; 
so Eq. 8 implies Eq. 9 but not conversely. They are only equivalent if the first equation in Eq. 8 is an equality. 
(When equality holds in Eqs. 8 and 9, the property is known as Strict Passivity.) 

We will say a system ( u , y) is Almost Strictly Dissipative (ASD) when there is some output feedback, 

U = G*y + U r , that results in the following being strictly dissipative: 

x = A c (x) + B(x)u r 

< A c (x) = A(x) + B(x)G*C(x) (10) 

J = C(x) 

Now if each component is ASD, then we have 


VI 'M (.V, ) + B, O, )G,C : <x : )] < -5, (x , ) + Z s F v A l (x, ,Uj) 

< . 7=1 (11) 

VV i B i (x i ) = Cf(x i ) 

where VV^- = gradient (V^ . Due to the interconnection terms, Eq. 11 is not necessarily Strictly Dissipative. 

However, in some circumstances, the interconnection terms have a special form and ASD is inherited when the 
system evolves. 

Suppose we have a pair of subsystems of the form: 

x t = A (x ) + sB i (x )u i + B* (x )IL 

<y i =C i (x i ) (12) 

X=C(A) 
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/ 

where i = 1, 2 and both subsystems 

U [ 

A 

A 

> 4 

\ / 

and 

u 2 

A '■ 

a 

> A 

\ 

have storage functions V- . We have the 

V 

u x 

La J, 

) V 

l U 2 J 

Lv 2 j 

) 


following result: 

Theorem 1: If the subsystems (u x , y x ) and ( U^y 2 ) are ASD and 


VV i B i (x i ) = Cf(x i )-i = \, 2 (13) 


then the resulting feedback connection, y l = U 2 and 

almost strictly passive. 


f 

A 


r at\ 


u x 


a 

III 

s 

A 

III 


V 

U 2 


La JJ 


U i=—y 2 9 will leave the composite system 


The physical connection of 
components in evolving systems is 
equivalent to the feedback connection of 
the admittance of one component to the 
impedance of the other component 3,5 . 
See fig. 2 for a schematic of the 
admittance-impedance connection 

between two evolving system 
components. The A superscript on the 
input and output variables is used to 
indicate that the component is in 
admittance-impedance form. For a 
nonlinear evolving system with 
components (« , y ) and (u^y^) 

connected in admittance-impedance 
form, Theorem 1 shows that ASD is an 
inherited property. 

Fig. 2. Admittance-Impedance feedback connection of evolving systems. 


Proof: See appendix. 


© 


u x - —u 2 

► 


r = 4 ( Vi I + sB ( Yi )//. 4- B A ( Vi ^\uf m 


AT “ C\(*i) 

I ' /*\ IVV y\ /*v T \ m m 


yf = c l A (x i ) component 1 


x, = v4,(x,) + sB^ixAu-) + 5 4 (x, )u A 

u 2 - y x 

A = c i( x i) 

y? = cpx 2 ) Component 2 




IV. Adaptive Key Component Control with Persistent Disturbance Mitigation 

This section describes the adaptive key component control with persistent disturbance mitigation. The key 
component is chosen to be component # 1 and will be modeled by the following nonlinear evolving system 
component with an external persistent disturbance: 

x { = \ (jc, ) + sB x (x { )u { + B x (x { )u x + T 1 (x l )u D 

< y, = C,(x,) (14) 

$ = (?(*) 

which will be respresented by the triple (A t (x, ), B l (x, ), Z? 4 (Xj ), Fj (x, ), C 4 (x t )) . All vector fields in Eq. 14 will 
have the appropriate compatible dimensions and be smooth in their arguments with a single equilibrium point at 0 in 
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a neighborhood U. The persistent disturbance input vector U D (V) is N D -dimensional and will be thought to come 
from the following Disturbance Generator: 


U D =&Z D 


z D =F z D ;z D (0) = z { 


(15) 


where the disturbance state Z D (f) is N D -dimensional. Such descriptions of persistent disturbances were first used 
to describe signals of known form but unknown amplitude 10 . For example, step disturbances yield 0 = 1 and F=0 

o r 


while sinusoidal disturbances can be described by 


© = [i o] 


and 


F = 


-a> D 0 


where the frequency CO D 


is known but the amplitudes are not known. 

Assume that the disturbance generator parameter F is known. In many cases this is not a severe restriction, e.g. 
step disturbances and sinusoidal disturbances are common forms of persistent disturbances. It turns out that it is 
better to rewrite the above in the following equivalent form: 


\ U D = 0 Zd 

Uz> = L(p D 


(16) 


where (p D is a vector composed of the known basis functions for the solutions of Eq. 15 and (L, 0) need not be 
known. The following example illustrates the basis function form for the disturbance generator: 

z D (t) = e Ft z D (0 ) 


= \<Pi it), cp 2 (0, ■ ■ .,<P„ D (o] z D (0) 

N D 

= Yj 7: r/PM)= L( PD 


i = 1 


Note that L is directly related to F via its columns but not to 0 . Some rearrangement of the entries in the 
columns of F is needed to create (p D . A simple example of the above is given by the following: 



'0 

f 


'o' 


T 

X = 

0 

0 

X + 

l 

U + 

l 


u D = a sin {co D t + b ) 

= a x sin (co D t) + a 2 cos (a> D t) 


(17) 


Assume component #1 


r u x 


r 

yi 

A 

9 

A 

vL m i 


L^i \) 


is ASD. Also let the following matching condition be true: 




( 18 ) 
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The matching condition means that there exists H* such that = ^(x^) . Let component #2 represent 

the remainder of the evolving system and assume it is strictly dissipative by choice of local controllers, that is: 


f x 2 = A 2 (x 2 ) + sB 2 (x 2 )u 2 

t>2 = C 2 (x 2 ) (19) 

The components are in admittance -impedance form so when they are joined U x = ~y 2 and U 2 = y l . The 
Adaptive Key Component Controller with Disturbance Mitigation works through the control input-output ports 
(u x , y x ) of component #1 and is given by: 

u t = G e y'\ + G n ( Pi) 

<G e =-yAyt) T r e ^e> o ( 20 ) 

Gp = — Ji (<p D ) Yd’Yd >0 

V 


This produces 



theorem shows: 


t-> 00 


■> 


0 with bounded adaptive gains ( G e ,G D ) as the following convergence 


Theorem 2: Assume that V x and V 2 are positive Vx ^ 0 and radially unbounded, and (^A(x),B(x),C(x)} are 

continuous functions of x and S ( X ) , above, is positive Vx ^ 0 and has continuous partial derivatives in x . 
Furthermore, assume: 


1) 


The conditions of Theo.l are satisfied; so that 


(Y 1 
u x 


Ml 

, A 


A 



L^i J ) 


is Almost Strictly Dissipative (ASD) 


2) The following matching condition is met: R(T 1 (x x )) Cl R(B X (x, )) 

3) (p D is bounded (or equivalently, F has only simple imaginary poles and no right half-plane poles) 


Then the adaptive Controller given by Eq. 20 produces X = 


Xo 


->0 with bounded adaptive gains 


( G e ,G D ) when Component 1 is joined with Component 2 into an evolved system and the outputs are 


y i =c i (x i > 


->o 


Proof: See Appendix. 

It should be noted that the above results might only hold on a neighborhood (0, r ( ) = jx. / ||x. || < j . Then 

the stability in Theo. 2 is only locally asymptotic to the origin. 

In the next section, an illustrative example will be given. 
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V. Illustrative Example 

Example 1 is a two component linear flexible structure evolving system, see fig. 3. The components of Example 
1 are stable when they are unconnected components, but the evolving system fails to inherit the stability of the 
components. This example will be used to demonstrate the inheritance and lack of inheritance of almost strict 

dissipativity in evolving systems. 

Component 1 Component 2 


V l=?l “ V 2 = <h V 3=«3 



Fig. 3. A two component flexible structure Evolving System. 


The dynamical equations for the 
components of Example 1 are: 

\m t q t =u x - s n k n (q l -q 2 ) 

=k>4i] T 

tn 2 q 2 — u 2 — s 12 k l2 ( q 2 — q 3 ) 
-k 22 (q 2 ~q 3 ) (21) 

< m 3 q 3 = u 3 - k 22 (q 3 —q 2 ) 

y 2 =[^ 2 ^ 2 ] t 

^3 = [^3 ’ ^3 ] 


with m l = 30 , m 2 = 1 , m 3 = 1 , k l2 = 4 , and k 22 = 1 . Example 1 has the following controllers: 


i x — -(Q.9^ + 0.1)^ 


ii 2 - — 


0.1 




f 0.2 s + 0.5 q 2 

V s J 


= -(0.6v + l)g 3 


( 22 ) 


When two components join to form an Evolved System, at their point of contact, their velocities are equal and the 
forces exerted are equal and opposite. If the two components are given by ( f x , V x ) and (/ 2 ,V 2 ), then the contact 
dynamics of the Evolved System can be represented by: 


j/i = -h 

(v, =q l =v 2 =q 2 


(23) 


This connection is modeled as the admittance of one component connected in feedback with the impedance of the 
other component, as described in section II. When we use this idea of the joining of two components of an Evolving 
System as the feedback connection of their admittance and impedance, we can apply Theorem 1 from above to 
determine whether almost strict dissipativity is inherited by the Evolved System. The subsystem components from 

Example 1 are stable in closed-loop form when they are unconnected, i.e., S X2 = 0 . When £ n = 1 , the system is 
fully evolved and it has a closed-loop eigenvalue at 0.17, resulting in an unstable Evolved System. 
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10 


Mass Displacement with Adaptive Control on Mass 1 


8 

6 

4 


— Mass 1 

Mass 2 

Mass 3 


A Simulink model was created to implement an adaptive 
key component controller for Example 1 as described in 
the previous section. Simulations were run in which the 



-4 - 
-6 - 

-8 1 1 1 1 1 1 1 1 1 

0 50 100 150 200 250 300 350 400 450 500 


connection parameter, £ n , ranged from 0 to 1, allowing 

the system to go from unconnected components to a fully 
Evolved System. The key component controller was able 
to maintain system stability during the entire evolution 
process when it used the input-output ports on mass 1 of 
component 1, see fig. 4. When component 1 was the key 

component, is ASPR. 


Time (seconds) 



Time (seconds) 

Fm. 5. Adantive kev coirroonent on mass 3. 


When the key component controller was located on 
component 2 and used the input-output ports on mass 3, 
stability was not maintained, see Figure 5. The adaptive 
key component controller was not able to restore stability 
on mass 3 because that system was not ASPR, i.e., it had 
nonminimum phase zeros at 0.005 15±0.2009i. 

A second study was performed where a proportional 
integral (PI) controller is used to control a rigid body. The 
plant has additional flexible body modes. The perturbation 
parameter 8 is used to connect or disconnect the flexible 
modes. In this way, a system with flexible modes that 
vary depending on the operating condition can be 
modeled. For the example described here, the system 
needs the adaptive key component control to maintain 
stability of the rigid and flex modes. Finally, one flex 
mode is defined that causes a non-minimum phase zero to 
appear in the plant when the feedback control is used. A 
residual mode filter is designed for this mode and 
successfully applied to remove the system that is causing 
the dissipativity condition to not be satisfied. The final 


paper will provide details of these models and results of the simulation. 


VI. Conclusion 

We have presented a result (Theorem 1) describing when an evolving system will inherit the almost strict 
dissipativity traits of its subsystem components. An example was given of successful inheritance of almost strict 
dissipativity and failed inheritance of almost strict dissipativity. This result allows a control system designer to 
determine a sufficient condition for an evolving system to use an adaptive key component controller to restore 
stability. We also presented a convergence result (Theorem 2) for an adaptive key component controller to restore 
stability during evolution and mitigate persistent disturbances. An illustrative example of a system with a rigid mode 
and flexible modes is used to demonstrate these ideas. 


Appendix 

Proof of Theorem 1: 

Fet (uf 9 yf ) be ASD. From (9) and (11), there exists G t such that 


VV4 C (*,) - vy [A i (x i )+ B^x^Cfix,)] 

< -5,. (x,. ) + SyV VfAq (x,. , u, , uf ) (40) 

vv i b; x (x) = c; x (x i ) t 
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If we connect (u t , V, ) in feedback with ( u 2 ,y 2 ) , then V, = u 2 and Uj = — y 2 and, use (12) and (13), then we 
have W 1 y4 12 (x l ,U l , ) = W l B l (Xj )ll l = C[ (Xj ) [ — >’ 2 J = — y { ,) ; 2 and similarly, 

^^ 2^21 ( A 2 , 1 / 2 , ^2 ) X * 


Let X - 



and, from (12), 


x = A(x) + B{x)u 


A (-^i ) "f" ^12^12 (^2 ) 

A 2 (X 2 ) + ^21^21 (^l ) 


+ 


A(x,) 

0 

u j 4 

0 

B 2 ( x 2 )_ 



y : 


yi 

y* 


= C(x) = 


Ci (Xj ) 

C 2 a (x 2 )_ 


(41) 


VVA(jc) = [Wj W 2 


with V = V l + V 2 , using (13) and £ . ( . = £ tj = £ from (3), 

A (X| ) + fi4 l2 (x 2 ) 

A 2 (x 2 ) + s 4 21 (X[) 

= W 1 A 1 (x 1 ) + f (-)A 2 ) + W 2 A 2 (x 2 ) + £ (y T 2 yj 
< -[S,(x,) + S 2 (x 2 )] + £{-y\ y 2 ) + f(j 2 
= ~S(x) 


and 


Vra(x) = [VF 1 VC 2 


C(V.) 

Cj{x 2 \ 


Bf(x 1 ) 


^(x 2 )J 


= C r (x) 


Therefore 
desired. End of proof. 


A 

r ai 


r*T\ 




y i 


A 

’}’a = 


V 

U 2 


l ^ 2 jj 


is ASD with output feedback 


~g; 

0 " 

yt 

+ 

Ar 

u x 

_ 0 

G 2_ 


1 

^ <N 

s 

1 


Proof of Theorem 2: 

Since the physical connection of Component 1 to Component 2 is equivalent to the feedback connection 

~ ~y 2 and U 2 ~ ^2 • 

/A A \ 

By Theo.l we have that the closed-loop system ( u x , y x ) below is ASD: 
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Rewrite (20), to obtain 


i, = A\(X\)~ sB , (x, )C 2 (x 2 ) + (jtj )(/,' 

! i 2 = A 2 (x 2 ) + £B 2 (x 2 )C l (x l )',0<£ <1 

$=<?(*) 


( 42 ) 


= G e yf + G D cp D = + G> 0 + AG// ; 


AG = G-G* =[AG e AG d ];jj = 
AG = G = -y?(yt) T r,y = 


y i 

^Pd 


Ye 0 

o Yd. 


>0 


Combining (21) and (22) yields: 

x, = A x c (x x ) — eB 1 (x 1 )C 2 (x 2 ) 

+{B*{x l )G' D +T l {x Y )OL\<p D + £ 1 a (x 1 )w 
with vv = AG/7 

< and G* = —H*6L from(19) 

and A c (X, ) - A (X, ) + B:' (X,. )G/C/' (X, ) 
x 2 = (x 2 ) + sB 2 (x 2 )C l (x l ); 0 < £ < 1 

>’, 4 =C,"(x,) 


(43) 


(44) 


Let V — V 1 + V 2 and we have: 


v=-s( X )+( y yw) 


Form V c 


1 

2 


tr(AGy 1 A G T ) and obtain from (22): 


V G =tr(AG/~ l AG T ) 
= -tr(yf(y*) T AG T ) 
= -tr(yf{ w) T ) 

= -{y i’ w ) 


Define: V (x, AG) = V (x) + V G (AG) and, from (24) and (25), we have: 


(45) 


( 46 ) 
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( 47 ) 


V(x,AG) = V(x) + V g (AG ) 

= -S(x) + (;y A ,w)-(;y A ,w) 
= -S (x) < 0 


This guarantees that all trajectories (x, AG) are bounded. If V (x,A G) is uniformly continuous or V (x, AG) is 
bounded, then Barbalat’s Lemma [10] yields: 


S(x) — >0 

v 7 t-> 00 


and the positivity and continuity of S (x) imply that X = 
Consider 


X 


V(x, AG) = -S(x) 

<|5(x)| 


dS(x) . 
x 


< 


< 


< 


dx 

8S(x) 


dx 

dS(x) 


dx 

dS(x) 


dx 


||A(x)|| + ||fl(x)||||w A 
| A(x)|| + ||fl(x)|| | AG|| ||C A (Xj)| 


-> 0 . 


which is bounded because (x,AG) is bounded, S(x) has continuous partial derivatives and 
( A(x), B(x\ C(x)) are continuous, and a continuous function of bounded x W is also bounded in t. 

So, y t = C.(x.) — — >0 because C-(x.) is continuous. End of Proof. 
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